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Abstract 



We consider certain boundary conditions supporting soliton solutions in the gene- 
ralized non- linear Schrodinger equation (AKNS). Using the dressing transformation 
(DT) method and the related tau functions we study the AKNSr system for the van- 
ishing, (constant) non-vanishing and the mixed boundary conditions, and their asso- 
ciated bright, dark and bright-dark N-soliton solutions, respectively. Moreover, we 
introduce a modified DT related to the dressing group in order to consider the free 
field boundary condition and derive generalized N-dark-dark solitons. We have shown 
that two— dark— dark— soliton bound states exist in the AKNS2 system, and three— 
and higher— dark— dark— soliton bound states can not exist. As a reduced submodel 
of the AKNSr- system we study the properties of the focusing, defocusing and mixed 
focusing-defocusing versions of the so-called coupled non-linear Schrodinger equation 
(r— CNLS), which has recently been considered in many physical applications. The 
properties and calculations of some matrix elements using level one vertex operators 
are outlined. 



Dedicated to the memory of S. S. Costa. 



1 Introduction 



Many soliton equations in 1 + 1 dimensions have integrable multi-component generalizations 
or more generally integrable matrix generalizations. It is well-known that certain coupled 
mTilti-ficld generalizations of the non-linear Schrodinger equation (CNLS) are integrable and 
possess soliton type solutions with rich physical properties (see e.g. [1, 2, 3, 4]). The model 
defined by two coupled NLS systems was earher studied by Manakov [5] . Another remarkable 
example of a multi-field generalization of an integrable model is the so-called generalized 
sine-Gordon model which is integrable in some regions of its parameter space and it has 
many physical applications [6, 7]. The type of coupled NLS equations find applications in 
diverse areas of physics such as non-linear optics, optical communication, biophysics, multi- 
component Bose-Einsten condensate, etc (see e.g. [1, 2, 3, 8, 9]). The multi-soliton solutions 
of these systems have recently been considered using a variety of methods depending on the 
initial-boundary values imposed on the solution. For example, the direct methods, mostly 
the Hirota method, have been applied in [10, 3] and in ref. [11] a wide class of NLS models 
have been studied in the framework of the Darboux-dressing transformation. Recently, some 
properties have been investigated, such as the appearance of stationary bound states formed 
by two dark-dark solitons in the mixed-nonlinearity case (focusing and defocusing) of the 
2— CNLS systems [4]. The inverse scattering transform method (1ST) for the defocusing 
CNLS model with nonvanishing boundary conditions (NVBC) has been an open problem for 
over 30 years. The two-component case was solved in [12] and the multi-component model 
has very recently presented in [13]. The inverse scattering [12, 13] and Hirota's method 
[10] results on N-dark-dark solitons in the defocusing 2-CNLS model have presented only 
the degenerate case, i.e. the multi-solitons of the both components are proportional to each 
other and therefore are reducible to the dark solitons of the scalar NLS model. 

The r— CNLS model is related to the AKNS^. system, which is a model with 2r depen- 
dent variables. As we will show below this system reduces to the r— CNLS model under a 
particular reality condition. Moreover, general multi-dark-dark soliton solutions of generali- 
zed AKNS type systems available in the literature, to our knowledge, are scarce. Recently, 
there appeared some reports on the general non-degenerate N-dark-dark solitons in the 
r— CNLS model [4, 14]. In [4] the model with mixed nonlinearity is studies in the context 
of the KP-hierarchy reduction approach. The mixed focusing and defocusing nonlinearity 
2-CNLS model presents a two dark-dark soliton stationary bound state [4] . Here we show 
that similar phenomenon is present in the AKNS2 system. However, three or more dark- 
dark solitons can not form bound states. In [14] general multi-dark-dark solitons have been 
obtained in the framework of the algebro-geometric approach of the CNLS^ model. 

In this paper we will consider the general constant (vanishing, nonvanishing and mixed 
vanishing- nonvanishing) boundary value problem for the AKNS^ (r = 1, 2) model and show 
that its particular complexified and reduced version incorporates the focusing and defocusing 
scalar NLS system for r = 1, and the focusing, defocusing and mixed nonlinearity versions 
of the 2-CNLS system, i.e. the Manakov model, for r = 2, respectively. We will consider 
the dressing transformation (DT) method to solve integrable nonlinear equations, which is 
based on the Lax pair formulation of the system. In this approach the so-called integrable 
highest weight representation of the underlying affine Lie algebra is essential to find soliton 
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type solutions (see [15] and references therein). According to the approach of [15] a common 
feature of integrable hierarchies presenting sohton solutions is the existence of some special 
"vacuum solutions" such that the Lax operators evaluated on them lie in some Abelian 
(up to the central term) subalgebra of the associated Kac-Moody algebra. The boundary 
conditions imposed to the system of equations must be related to the relevant vacuum 
connections lying in certain Abelian sub-algebra. The soliton type solutions are constructed 
out of those "vacuum solutions" through the so called soliton specialization in the context of 
the DT. These developments lead to a quite general definition of tau functions associated to 
the hierarchies, in terms of the so called "integrable highest weight representations" of the 
relevant Kac-Moody algebra. However, the free field boundary condition does not allow the 
vacuum connections to lie in an Abelian sub-algebra; so, we will introduce a modified DT 
in order to deal with this case. We consider a modified DT relying on the dressing group 
composition law of two successive DT's [29] . The corresponding tau functions will provide 
generalized dark-dark soliton solutions possessing additional parameters as compared to the 
ones obtained with constant boundary conditions. 

Here we adopt a hybrid of the DT and Hirota methods [7, 18] to obtain soliton solutions 
of the AKNS system. As pointed out in [17], the Hirota method presents some drawbacks 
in order to derive a general type of soliton solutions, noticeably it is not appropriate to 
handle vector NLS type models in a general form since the process relies on a insightful 
guess of the functional forms of each component, whereas in the group theoretical approach 
adopted here the dependence of each component on the generalized tau functions become 
dictated by the DT and the solutions in the orbit of certain vacuum solutions are constructed 
in a uniform way. The generalized tau-functions for the nonlinear systems are defined as 
an alternative set of variables corresponding to certain matrix elements evaluated in the 
integrable highest-weight representations of the underlying affine Kac-Moody algebra. In this 
way we overcome two difficulties of the methods. The Hirota method needs the tau-functions 
and an expansion for them, but it does not provide a recipe how to construct them. That is 
accurately solved through the DT method, which in its turn needs the evaluation of certain 
matrix elements in the vertex operator representation. We may avoid the cumbersome matrix 
element calculations through the Hirota expansion method. Since this method is recursive 
it allows a simple implementation on a computer program for algebraic manipulation like 
MATHEMATICA. 

The paper is organized as follows. In section 2 we present the generalized non-linear 
Schrodinger equation (AKNS2) associated to the homogeneous gradation of the Kac-Moody 
algebra s/3. The 2— component coupled non-linear Schrodinger equation (CNLS2) is defined 
as a particular reduction by imposing certain conditions on the AKNS2 fields. In section 
3 we review the theory of the DT, describe the various constant boundary conditions and 
define the tau-functions. In 3.1 we present the modified DT associated to the dressing group 
suitable for fee field boundary conditions. In section 4 we apply the DT method to the 
vanishing boundary conditions and derive the bright solitons. In section 5 we consider the 
non-vanishing boundary conditions and derive the AKNS^ (r = 1, 2) dark solitons. Moreover, 
the general N-dark-dark solitons with free field NVBC of the AKNS2 arc derived through 
the modified DT method. In 5.3 we discuss the dark-dark soliton bound states. In 6 the 
dark-bright solitons are derived associated to the mixed boundary conditions. Finally, we 
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have included the sl{2) and s/(3) afiine Kac-Moody algebra properties, in appendices A and 
B, respectively. In C the matrix elements in s/(3) have been computed using vertex operator 
representations. 



2 The model 

The AKNS2 model can be constructed in the framework of the Zakharov-Shabat formalism 
using the hermitian symmetric spaces [19]. In [20] an affine Lie algebraic formulation based 
on the loop algebra g(E)C[\, A^^] of (7 G s/(3) has been presented. Here instead we construct 
the model associated to the full affine Kac-Moody algebra Q = sl{3) with homogeneous 
gradation and a semi-simple element . So, the connections are given by 



A = Em + Y^nE'£' + T.^7E'y,^ + <PiC, (2.1) 

i=l i=l 

B = Ei^^+j:nE^;^+j:KE^2i+±d^nEfy±d.^rE^^l- 

1=1 1=1 1=1 1=1 

[v^+v^r - ^i]{Hi'^ + ) - [^^^2 - ^2]hP - 

^t^2Efl0, - ^2+^r4?-/3. + <I>2C, (2.2) 

where \E'^, (pi and 02 ai'e the fields of the model. Notice that the auxiliary fields 0i and 02 
lie in the direction of the affine Lie algebra central term C. The existence of this term plays 
an important role in the theory of the so-called integrable highest weight representations of 
the Kac-Moody algebra and they will be important below in order to find the soliton type 
solutions of the model. The connections (2.1)- (2.2) conveniently incorporate the terms with 
the constant parameters i^i 2 in order to take into account the various boundary conditions, 
so, these differ slightly from the ones in [21, 20]. 
In the basis considered in this paper one has 



E^'^-l^aHi^^ (2.3) 



a=l 



The /3i in (2.1)-(2.2) are positive roots defined in (B.24). Notice that the connections fie in 
the subspaces 

Aego + gi, Be 90 + 91+92. (2.4) 

The zero curvature condition [dt —B ,8^ — A] — provides the following system of 
equations 



= +aX+_2[^^r+^T-l(/3,.d)]*+, i = l,2 (2.5) 
dt^T = +2[E*.^*7-^(^-^)]*^^' (2-6) 
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dt4>i-dj2 = 0; (2.7) 

2 

Q = ^aA, ( A. j1) = 2^1 + ^2, (^.jl) = 2^2 + 

i=l 

Notice that the auxihary fields 0i,2 completely decouple from the AKNS2 fields ^f. The 
integrability of the system of equations (2.5)-(2.6), for Vti^2 = 0, and its multi-Hamiltonian 
structure have been established [21]. The system of equations (2.5)-(2.6) supplied with a 
convenient complexification can be related to some versions of the so-called coupled non- 
hnear Schrodinger equation (CNLS) [5, 1, 2, 3]. For example by making 

t^-it, = -n Gi *7 = (ji ^i, (2.8) 

where ★ means complex conjugation, // e IR+, (jj = ±1, we may reduce the system (2.5)-(2.6) 
to the well known 2— coupled non-hnear Schrodinger system (2 — CNLS) 

idtiPk + dliPk + ^fi (^a^\iPj\^-^\{f3k.fl)\^iPk^0, k^l,2. (2.9) 

The parameter fi > represents the strength of nonlinearity and the coefficients aj define 
the sign of the nonlinearity. The system (2.9) can be classified into three classes depending 
on the signs of the nonhnearity coefficients (Tj. For (Ji = (T2 = 1 this system is the focusing 
Manakov model which supports bright- bright solitons [5]. For (Ji = (T2 = — 1, it is the 
defocusing Manakov model which supports bright-dark and dark-dark solitons [10, 12, 22]. In 
the cases ai = — (T2 = ±1 one has the mixed focusing-defocusing nonlinearities. In this case, 
these equations support bright-bright solitons [23, 3], bright-dark solitons [24]. This mixed 
nonlinearity has recently been considered in [4] through the KP-hierarchy reduction method 
and dark-dark sohtons have been obtained. The system (2.9) has also been considered in 
the study of oscillations and interactions of dark and dark-bright solitons in Bose-Einstein 
condensates [25, 26]. The multi-dark-dark solitons in the mixed non-linearity case are useful 
for many physical applications such as nonlinear optics, water waves and Bose-Einstein 
condensates, where the generally coupled NLS equations often appear. 

The focusing CNLS system possesses a remarkable type of soliton solution undergoing a 
shape changing (inelastic) collision property due to intensity redistribution among its modes. 
In this context, it has been found a novel class of solutions called partially coherent solitons 
(PCS) which are of substantially variable shape, such that under collisions the profiles remain 
stationary [1, 27, 28]. Interestingly, the PCSs, namely, 2-PCS, 3-PCS, r-PCS, are special 
cases of the well known 2-soliton, 3-soliton,..., r-soliton solutions of the 2-CNLS, 3-CNLS,..., 
r-CNLS equations, respectively [2, 3]. So, the understanding of the variable shape collisions 
and many other properties of these partially coherent solitons can be studied by providing 
the higher-order soliton solutions of the r-CNLS (r > 2) system considered as submodel 
of the relevant AKNS^. We believe that the group theoretical point of view of finding the 
analytical results for the general case of A'"-sohton interactions would facilitate the study 
of their properties; for example, the asymptotic behavior of trains of solitonlike pulses 
with approximately equal amplitudes and velocities, as studied in [16]. Notice that the set 
of solutions of the AKNS model (2.5)-(2.6) is much larger than the solutions of the CNLS 
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system (2.9), since only the solutions of the former which satisfy the constraints (2.8) will 
be solutions of the CNLS model (2.9). This fact will be seen below in many instances; so, 
we believe that the AKNS soliton properties with relevant boundary conditions deserve a 
further study. 

We will show that the three classes mentioned above, i.e. the focusing, defocusing and 
the mixed focusing- defocusing CNLS model can be related respectively to the vanishing, 
nonvanishing and mixed vanishing-nonvanishing boundary conditions of the AKNS model 
in the framework of the DT approach adapted conveniently to each case. The (constant) 
non- vanishing boundary conditions require the extension of the DT method to incorporate 
non-zero constant vacuum solutions. Therefore, the vertex operators corresponding to the 
vanishing boundary case undergo a generalization in such a way that the nilpotency property, 
which is necessary in obtaining soliton solutions, should be maintained. Using the modified 
vertex operators we construct multi-soliton solutions in the cases of (constant) non-vanishing 
and mixed vanishing-nonvanishing boundary conditions. The free field NVBC requires a 
modification of the usual DT of [15] by considering two successive DT's in the context of the 
dressing group [29]. However, the same vertex operator generating the dark-dark solitons in 
the (constant) NVBC will be used in the free field NVBC case with a modified tau functions. 



3 Dressing transformations for AKNS2 

In this section we summarize the so-called DT procedure to find soliton solutions, which 
works by choosing a vacuum solution and then mapping it into a non-trivial solution, follow- 
ing the approach of [15]. For simplicity, we concentrate on a version of the DT suitable for 
vanishing, (constant) non-vanishing and mixed boundary conditions of the AKNS2 model 
(2.5)-(2.7). The free field boundary condition requires a modified DT which is developed in 
subsection 3.1. So, let us consider 

lim^^^-p^; 012^-0; = const. (3.1) 

We may identify the NVBC (3.1) to certain classes of trivial vacuum solutions of the 
system (2.5)-(2.6): 

1) the trivial zero vacuum solution 

*,>ac = Pf = 0, J = 1,2 (3.2) 

2) the trivial constant vacuum solution 

*,>ac = Pf7^0, J = 1,2. (3.3) 

Notice that (3.3) is a trivial constant vacuum solution of (2.5)-(2.6) provided that the 
expression [I]|=i Pj p] — |(/9i.n)] vanishes for any i = 1, 2; which is achieved if = O2 = 
implying p^ p^ = ffi. 

3) the mixed constant-zero (zero-constant) vacuum solutions 

0^t.ac = P^^O, vl/±_=p±=o, (3.4) 
^^)^l,vac = pf = Q, ^tvac = pt^Q- (3-5) 
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The first mixed trivial solution (3.4) requires 2piPi = 2Qi + Q2, whereas for the second 
trivial solution (3.5) it must be 2p2P2 = 21^2 + l^i. 

The connections (2.1)-(2.2) for the above boundary conditions (3.1) take the form 

^.ac ^ E(^^+ptEf^+p^E^\+ptE^;^+p^E% (3.6) 
Bvac ^ E^^)+ptEf^+p-,E%^+ptEf^+p2E%-ptp2Efl,^-ptp-,Efl,^- 

(ptpT - n,){H?^ + i/f ) - iptp2 - n2)Hi'\ (3.7) 

Notice that ^A"^'^ , B'""''^] — 0. These connections are related to the group element "^^^^ 
through 



where 



^(0) _ gXA'""=+tB'"" 



(3.8) 



(3.9) 



The DT is implemented through two gauge transformations generated by Q± such that 
the nontrivial gauge connections in the vacuum orbit become [15] 



(3.10) 
(3.11) 



where 



= exp ( ^ a_ J , QX = M-^N; M = exp ((7„) , TV = exp ( ^ , (3.12) 

\n>0 / \n>0 / 



where [I?, (7„] —ncTn- Therefore, from the relationships 
one has 



(3.13) 



(3.14) 



where h is some constant group element. 

One can relate the fields \E'f , 0i and 02 to some of the components in cr„. One has 



A — (7_i,£^^^^ + terms of negative grade. 



(3.15) 



+ 2 



+ 



(1) 



i=l 



+ 



j=i 



+ terms of negative grade 



(3.16) 
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Taking into account the grading structure of the connection A in (2.4) we may write the 
cr^s in terms of the fields of the model. In order to match the zero grade terms of the both 
sides of the equation (3.15) one must have 



(3.17) 



i=l 



a=l 



In the equation above, the explicit form of cr" ^ in terms of the fields "^f can be obtained by 
setting the sum of the (—1) grade terms to zero. Nevertheless, the form of the a°ii will not 
be necessary for our purposes. 

Following the above procedure to match the gradations on both sides of eqs. (3.15)-(3.16) 
one notices that the a-n 's with n > 1 are used to cancel out the undesired components on 
the r.h.s. of the equations. 

From the equations (3.12)-(3.13) and (3.14) one has 



(A|M-i 



A 



^(0)/j,^(0)- 



(3.18) 



where (A| and A'^ are certain states annihilated by Q< and respectively. Defining 
= E + E ^o'-E^-l + E + r)C (3.19) 



a>0 



a>0 



a=l 



and choosing a specific matrix element one gets a space time dependence for the field r] 



(Aol |Ao) 

To, 



(3.20) 
(3.21) 



where we have defined the tau function tq. Next, we will write the fields "if in terms 
of certain matrix elements. These functions will be represented as matrix elements in an 
appropriate representation of the affine Lie algebra sl{3). We proceed by writing the eq. 
(3.14) in the form 



exp - E ^-n |Ao) = |A„) T,-\ 



n>0 



(3.22) 



where the eqs. (3.12), (3.19) and (3.20)-(3.21) have been used. 

Then the terms with grade (-1) in the both sides of (3.22) can be written as 



-cr-i \Xo) 



(-1) 



|A„ 



To (X, t) 



(3.23) 



or equivalently 



fE(*.^ - Pt)E^,f - E(*r - P7)E^ - E <^-iHtA I Ao) = 

\i=l i=l a=l ) 



^(0)/^^(0)-l 



(-1) 



|Ao) 



To {X, t) 



(3.24) 
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Acting on the left in eq. (3.24) by E^j^^ and taking the relevant matrix element with the 
dual highest weight state < Xo\ we may have 

'^t^pt + — and ^r = pr-!L; i^i,2. (3.25) 

To To 

where the tau functions t^, Tq are defined by 

r+(x, t) ^ (A„| E% |A„) , (3.26) 

Tr{x,t) ^ (A,|4;) |A„) , (3.27) 

ro(a;,t)^(A„|[*(°)/i*(oH^ IA„). (3.28) 

L J (o) 

Notice that in order to get the above relationships we have used the commutation rules 
for the corresponding affinc Kac- Moody algebra elements, as well as their properties acting 
on the highest weight state (sec Appendix B). 

According to the solitonic specialization in the context of the DT method the soliton 
solutions are determined by choosing suitable constant group elements h in the eqs. (3.26)- 
(3.28). In order to obtain n— soliton solutions the general prescription is to parameterize the 
orbit of the vacuum as a product of exponentials of eigenvectors of the operators ei {si = 
A"""", £2 = 5^"^) defined in (3.6)-(3.7); i.c h = IVl^ie^\ where [ei, Fi] = \\Fi, , such that 
(Fj)"* 7^ only for m < rrii, rrii being some positive integer. The relationships between 
DT, solitonic specialization and the Hirota method have been presented in [15] for any 
hierarchy of integrable models possessing a zero curvature representation in terms of an 
affine KacMoody algebra. The DT method provides a relationship between the fields of the 
model and the relevant tau functions, and it explains the truncation of the Hirota expansion. 
The Hirota method is a recursive method which can be implemented through a computer 
program for algebraic manipulation like MATHEMATICA. On the other hand, the DT 
method requires the computation of matrix elements as in eqs. (3.26)-(3.28) in the vertex 
operator representations of the affine KacMoody algebra. Actually, these matrix element 
calculations are very tedious in the case of higher soliton solutions. 

3.1 Free field boundary conditions and dressing group 

As a generalization of the constant NVBC (3.3) consider the free field NVBC 

hm ^ pf {x, t) ^ e"f * (3.29) 

where pf, af, hf are some constants such that of — —a~; bf — —h~. Notice that (3.29) 
is a free field solution of (2.5)-(2.6) such that hf = ±(a,^)^ T 2Aj, for Aj = [E5=iP/p~ - 

K^j.Q)] = const. However, a direct apphcation of the DT approach of [15] is not possible 
since the relevant connections A^"*^, B'""''^ do not belong to an abelian subalgebra (up to the 
central term). So, in order to consider the NVBC (3.29) we resort to the dressing group 
composition law of two successive DTs. In fact, the relevant connections can be related 
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to the A'""''^, B'""-'^ in (3.6)-(3.7) through certain gauge transformations generated by 6± such 
that 

yvac ^ eux,t)v;""'[eux,t)]-^ + d^ieux,t)[eux,t)]-^; Vi = V2 = (3.30) 

where i = 1, 2, xi = x, x-i = = A"""", V2 = -B^«^ and Q% e §L(3). One must have 

^(0) ^ ^(0) ^ e^(a;,t)^(°) = e^(a:,t)^(°)y, = (3.31) 

is an arbitrary constant group element, and 

Qfl ^ e^°eS">o^", e?. = e->^°eSn>o>^-", [D , ^n] = riXn- (3.32) 
Prom (3.30) one so 

Xo = X:<l/3.+Xo-^i%. + Ex:^f +XC' (3.33) 

a=l 

x-1 = -E[/^(^,^)-P.^]47'^ + E[r(^,i)-pr]^i",f + Ex^^r^^- (3.34) 

i=l i=l a=l 

Therefore, a modified DT procedure can be implemented with [29] 

= ^(o)/i'[^(o)]-i = el(x, i)*(°)/i[*(°)]-^[e^(x, t)]-^- h = gh', (3.35) 

through 

A' = ei'M('''^'=)[ei'^]-^ + d^el'iel']-^; b' = ei'^B('''^<^)[ei'^]-^ + dsi'[Qi']-\ (3.36) 

So, the equation (3.35) can be used instead of (3.14) in order to derive general dark-dark 
solitons with free field NVBC. 

4 Bright solitons and vanishing boundary conditions 

For simplicity we first apply the DT method to the VBC (3.2) and show the existence of 
bright- bright soliton solutions [30]. The connections A^"-'^ = E^^\ B'""-'^ = E^'^^ are related to 
the group element as 

^(^«-) = a,^o[*o]"'; = 9t^o[^o]"'; ^0 = e^^*''+*^*''. (4.1) 

In order to obtain soliton sohitions the simultaneous adjoint cigcnstates of the elements 
E^^\ E^"^^ play a central role. In the case at hand one has the eigcnstates 

+0O +0O 

Fj- E E Pk47^-^ j,k^h'2 (4.2) 

n=—oo n=— oo 

such that 

xE^^^ + tE^^\ Fj\ = -<Pj{x, t)Fj] ipj{x, t) = Uj {x + Ujt) (4.3) 

'xE^^^ +tE^^\ Gj] = r]j{x,t)Gj, rij{x,t) = pj {x + pjt) (4.4) 
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4.1 J component one bright-soliton solution 

Consider the product 

/i = e"^i-^-'ie''^2^^2, (4.5) 

where ji and j2 are some indexes chosen from {1,2}. Using the nilpotency properties of Fj^ 
and (see Appendix C) one gets 

'^oh%'] = (l + e-^^ia,,F,,)(l + e^^-^6,,G,J (4.6) 

= 1 + e-"^" Uj, Fj^ + e''« bj^ Gj^ + aj, bj^ e"'^^! e^^'^ Fj^ Gj^ , (4.7) 

with (pj^ and r^jj given in (4.3)-(4.4). The corresponding tau functions become 



Tn = 1 + a b C- • e~'^-'ie''-'2 C ■ = '^llfll A. . (as.) 

K^h - Ph) 

— 3i,j2^j2Pj2^^^^ 1 '^i ~ ^i,h^h^h^ '^n ^ (4-9) 

where the matrix element Cj^j^ has been presented in (C.2). In order to construct one-soliton 
solutions we must have ji — j2 = j in (4.8). Therefore one gets ^ 

^+ = -L^ qrr = — i = j- (4 10) 

' 1 + aibiCne-Vi e^i' ' 1 + 6, e-^» e''» ' ^ ^ 

*± = 0, i^j (4.11) 

Impose the relationships p* = — z/j, 6* = —/laaj i^j^aj G €; so from (4.8) one has Cjj = 
_fj^\2 rjj-^g equations (4.10)-(4.11) with the complexification (2.8) provide a solution of 



the CNLS system (2.9) 

f — M^e^^^^sech ((piR + ^) , i=j 

where = ^fe^' = " ^^jO = V^jR + ^'/^j^'^- 

The solution (4.12) for the j'th component is known as a 'bright soliton' in the context 
of the scalar non-linear Schrodinger equation (NLS). 



4.2 1-bright-bright soliton solution 

The main observation in the last construction of the j'*'*— component one-soliton is that it 
has been excited by the group element h in (4.5) such that ji = j2 = j- So, in order to excite 
the two components of {i = 1,2) and reproduce a bright soliton for each component let 
us consider the group element 

;^^gaiFig6iGiga2F2g62G2^ (4.13) 
^If ji ^ j2 in (4.8)-(4.9) one still has certain trivial solutions, since in this case Cj^j^ = implying tq = 1 
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where the exponential factors contain F's and G's of type (C.l). The tau functions become 



r+ = bjpje''^, rjj^ pj{x + pjt)+poj, j = l,2. (4.14) 
T~ = ajUje~'^% (fij = Uj {x + Ujt) + ugj, j = 1, 2. (4-15) 
To = l + ai6iCne-'^ie^i +a262C22e-'^^e^^ (4.16) 

where Cjj = -^^^^z^- Let us set b* = —fio-jaj, p* = —Uj = — z/i,then rj* = —ipj = —ipi = 
~{^iR + i^ii), ^ ^^) ill the relations (4.14)-(4.16). Therefore, the vector soliton 

solution of the 2-CNLS eq. (2.9) arises 



j uiRsech(^^,R + ^^ e^'^", (4.17) 



Notice that this 1-bright-bright soliton possesses 3 arbitrary complex parameters, namely, 
Oi, 02 and = I'm + ii'u- This solution in the case of mixed nonlinearity cxi = — (J2 = 1, may 
have a singular behavior if the sum {(Tj\ajVj\'^)^^'^ in the denominator of the expression of 
Ai vanishes, such that the soliton amplitude in (4.17) diverges. The A^— bright-bright soliton 
requires the generalization of the group element in (4.13) as h — e"i''^ie''^*^^...e"^^^e''^'^^ 



5 AKNS^ (r = 1, 2), NVBC and dark solitons 

In this section we tackle the problem of finding dark soliton type solutions of the system (2.5)- 
(2.6). The associated coupled NLS model (2.9) with nonvanishing boundary conditions has 
been considered in the framework of direct methods, such as the Hirota tau function approach 
(see e.g. [10, 31, 22, 32]), and recently in the inverse scattering transform approach [12, 13]. In 
the last approach the relevant Lax operators have remarkable differences and rather involved 
spectral properties as compared to their counterparts with vanishing boundary conditions 
(see [33] and references therein), e.g. in the NVBC case the spectral parameter requires the 
construction of certain Riemann sheets [33, 34, 35, 36]. So, it would be interesting to give 
the full Lie algebraic construction of the tau functions and soliton solutions for the system 
(2.5)-(2.6) with NVBC. 



5.1 AKNSi: N-dark solitons 

For simplicity, firstly wc describe the entire process for the system (2.5)-(2.6) with just two 
fields So, in order to study a NVBC for the system s/(2)-AKNS consider the Lax pair 

A = + + + (/)iC, (5.1) 

d:,^!^Ef - d^^~E^^^ - 2 - p+p-) + 02C, (5.2) 
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where ^~ are the fields of the model {p^ = constant). In this case one considers the 
generator E^^^ = H^''\ The 0i and 02 are introduced as auxiliary fields. Therefore the 
equations of motion suitable to treat NVBC become 

= - 2 (5.3) 

dt^- = -9^*- + 2(*+*- -p+p-)*-, (5.4) 
dt(t>i-d,(t>2 = 0. (5.5) 

We will show that the s/(2)-AKNS model, through a particular reduction, contains as 
sub-model the scalar defocusing NLS system 

iatV' + a^V'-2 (|v^|2-p2)v^ = o. (5.6) 

This equation is suitable for treating nonvanishing boundary conditions (NVBC) [36, 33, 



34] 



In fact, this boundary condition is manifestly compatible with the eq. of motion (5.6). 
The defocusing NLS equation (5.6) with the boundary condition (5.7) is exactly integrable by 
the inverse scattering technique [37] . This model has soliton solutions in the form of localized 
"dark" pulses created on a constant or stable continuous wave background solution. 

In the system of eqs. (5. 3) -(5. 4) consider the transformation 

t — >■ —it, X — >■ —ix, (5.8) 
^ ^*^e^^ (5.9) 

where the factor e^^ is introduced for later convenience. Furthermore, the identification 

V; = = (*-)*, (5.10) 

where the star stands for complex conjugation, such that ^'o^^'o ^ —p^ allows one to 
reproduce the defocusing NLS equation (5.6). 

Let us take as vacuum solution of (5.3)-(5.4) the constant background configuration 

'^^^p^ = p€^^, 01,2 = 0,, p,€ = const. (5.11) 

Therefore the gauge connections (5.1)- (5. 2) for the vacuum solution (5.11) become 

A""" = El = E^^^ + p+ Ef + p- E^^^ (5.12) 
B^'"' = £2 = ^^'^ + P+£^+ ^ + ^ (5.13) 

We follow eq. (3.8) in order to write the connections in the form 

^ivac) ^ a.^Wj^Wj-l, i^(-) = a,*il[*il]-\ (5.14) 



12 



with the group element 

= e^' 

= exp [xiE^'^ + p-^Ef + p-£;i°)) + + p+£;« + p-£;i'))] (5.15) 



Vl>(0) _ a;£i+te2 



Let us emphasize that this group element differs fundamentally from the one associated 
to the case with VBC. First, the difference originates from the constant boundary condition 
terms added to the relevant vacuum gauge connections. However, they must be related since 
lini — )■ ^^0, where is the group element in the VBC case (4.1). Second, the vacuum 

connections £i_2 in the DT procedure will require another eigenvectors under the adjoint 
actions in analogy to eqs. (4.3)-(4.4), as we will see below. 

Another proposal for this group element has been introduced in [38]. There it 

has been considered a group element inspired in the inverse scattering method, possessing 
double-valued functions of the spectral parameter A. This fact motivates the introduction of 
an afRne parameter to avoid constructing Riemann sheets. This construction involves some 
complications when used in the context of the DT method. The main difficulty arises in the 
computation of the matrix elements associated to the highest weight states in order to find 
the relevant tau functions. However, the group element given in (5.15) will turn out to be 
more suitable in the DT procedure. In fact, a similar group element has been proposed in 
[39] to tackle a NVBC soliton solutions in the so-called negative even grade mKdV hierarchy. 

The DT procedure in this case follows all the way verbatim as in section 3 adapted to the 
sl{2) case. It follows by relating the relevant connections (5.1)-(5.2) with the connections in 
eqs. (5.12)-(5.13) corresponding to the vacuum solution (5.11). 

Next, we will write the fields in terms of the relevant tau functions. According 
to the development in section 3 these functions will be written as certain matrix elements 
in an integrable highest weight representation of the affine Lie algebra sl{2). So, it is a 
straightforward task to get the following relationships 



where the tau functions r^, tq are given by 



+ 



r 

p+ — and ^~ 

To 



(1) 



^nvbc"'^ 



P 



(0)-l 
nvbc 



(-1) 



r 

To 

[A. 



^(0) ;,^(o)-i 



nvbc 



-1) 



|Ao), 



nvbc nvbc 



To = (A 

with the group element '^^^Ibc given by (5.15). 



|Ao), 



(5.16) 

(5.17) 
(5.18) 

(5.19) 



5.1.1 AKNSi and l-dairk soliton of defocusing scalair NLS 

Let us consider the group element 



, (?=1,2) 



(5.20) 
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where a'^ is some constant. 

The vertex operators V'^ and V'^ in (A. 10), respectively, satisfy 

[£1,^ = 27^^^ = (-l)«-' 27(72 -p2)i/2v^<z^ q^i^2. (5.21) 

Notice that these eigenstates exhibit a more complex structure in comparison to their 
counterparts Fj and Gj corresponding to the VBC as in (4.3)-(4.4). Therefore one has 

[xe^ + te2 , V^i^,, p)] = [27, {x + {-iy-\t)] F«(7,, p), (5.22) 
where Vq — — p^. Denoting ipq = 27^ [x + {—lY~^Vqt) one can write 



nvbc novo 



exp {e'^'^aqV^) (5.23) 
1 + e^-'aqV\ (5.24) 



where we have used {V^Y = 0) for n >2, that is, is nilpotent (A.15)-(A.16). 

So, in order to find the explicit tau functions in (5.17)-(5.19) it remains to compute 
the relevant matrix elements. Using the properties presented in the Appendix A and the 
eqs. (3.26)-(3.28) one gets the tau functions rj) = 1 + d'^^ e^"; rj^ = s±(«)p±e'^«, where 

the matrix elements c^") = a^^^ {Xo\ \\o) , = a^^^ {p^)'^ (Ao| E^^W'i have been 

considered. They are given in the Appendix A, eqs. (A.11)-(A.12), respectively. Therefore, 
the relations (5.16) provide 

*±('')(x,t) = p^[l ± ^] ± ^^^tanh{7,[a:+ (-l)^^V] + ^-}. (5.25) 

The solution in the case q = 1 allows the imposition of the conditions (5.10) in order 
to satisfy the defocusing NLS equation (5.6). So, let us make the transformation (5.8) and 
the complexification 71 — )■ ^7, p — )■ ip, ^ = \/ p"^ — 7^; (p > 7). Choose ai = 
—i^jv. Moreover, considering the transformation (5.9) one gets the functions satisfying 
the complexification condition (5.10). Therefore, the function i){x^t^ = — ie^^"*", where 

= ie~2{p+ ( jj^l^ti ) i+e-^c'^^t) }; ^{^-i t) = 2^{x — v^f) +7o, obtained in this way is a solution 
of the defocusing NLS equation (5.6). In fact, the function -0 can be written as 

il){x,t) = e^(°^+'^/2){^tanh[7(a; - v^t) + 70/2] + ^ - ipe"^"}, a = tan-^(^), (5.26) 

which is the dark soliton solution of the defocusing NLS equation (5.6) (see e.g. [31]). Notice 
that this solution satisfy the NVBC lim 1-0(2;, i)!^ ^ p^. 

Similarly, one can consider the case q — 2. In this case one gets = ie^^p 

^e='=**tanh(7a; + 'yv^t) + 1 ^ 2^^^^^ ' ^ = tan~^(v^/7). However, they do not satisfy the 
condition (5.10), hence these functions would not be solutions of (5.6). 

From the previous constructions it is clear that the two types (each one associated to 
the index q = 1,2) of N-dark solitons of the AKNSi model will be associated to the group 
element 

f^q ^ e<v<'(7,\p)ga2y<'(72,p)___ga^y«(7f ,p)^ g = 1, 2 (5.27) 
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5.2 AKNS2: N-dark-dark solitons 

Consider the non-vanishing boundary condition (NVBC) for the system (2.5)-(2.6) in the 
form (3.1) and its associated constant vacuum solution (3.3). The vacuum connections (3.6)- 
(3.7) associated to this trivial solution commute [^""'^' . j^^a-cj = if one takes Qi = Q2 ~ 
UpiPi + P2P2)- Consider the notation A^"^ = si , B"'"' = 82- 

l-d£irk-deirk soliton. Let us choose the group element /i* = ga'w«(fe)^ j^q i-,ging 

some parameters, in the eqs. (3.26)-(3.28), such that 

[ei,]^^] = k'^W'^; q=l,2; (5.28) 
[£2, W] = w'^W'i (5.29) 

where 



So, one has 



e,k^^{kir-AptpT-^ptP2, e, = {-iy-' (5.30) 



[xei + tEi^W^] = {k'^x + wH)W^ (5.31) 



The vertex operator (see (C.7)) is associated to the 1-dark-dark sohton and it can be 
shown to be nilpotent. So, using (3.26)-(3.28) the following tau functions correspond to the 
element h'^ given above 

ri«) = i + c(«)e^''-+-''*, (5.32) 
rf(«^ = p±s±We'='^+"''*, i = l,2; (5.33) 



where c^^) = < Xo\Wi\Xo >, s^^^^ = {pf)~^a'i < XolE^^^W^lXo >■ The above tau 
functions provide a 1-dark-dark solution for the fields of the system (2.5)-(2.6), such 
that the following relationships hold between the parameters 

^ ' ' g+ {q) — g- {q) ^ ' ^ ' 

Then, the relations (3.25) provide 



± 

^ 2 



Ml) - P. 



(1 + (|/'')±i) + (-1 + (1/^)^^) tanh{(A;% + wH + /o^ c^-^Vs}] , (5.35) 



where j = 1, 2; = pr(^. The condition c*^^^ > requires that > 1 and < 0. 

Taking into account (C.7) one can compute the relevant matrix elements, which define 
s± (9) ^ in order to get 

(6.36) 

- e,v'(*')"-'tEjP>7 



Notice that for any value of the index g, if Ylij P^Pj > the last relationship implies 
y'' <0 and if pjpj < then y'^ > 0. 
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Some remarks are in order here. 

First, the sohtons associated to the pair of components {'^t ^ ^2) ^-^d (^'J" , ^'2), re- 
spectively, arc proportional, so they are degenerate and reducible to the single dark soliton 
of the sl{2) AKNSi model (5.25). Whereas, the solitons associated to the pair of compo- 
nents {"^t , ^7)' respectively, arc not proportional, so they are non-degenerate presenting 
in general different degrees of 'darkness' in each component. In this way the solution in 
(5.35) presents a non-degenerate 1-dark-dark soliton in the components, say {"iff , ) of 
the AKNS2 model. A nondegenerate 1-dark-dark sohton in the defocusing 2-CNLS has been 
given in [22, 10] through the Hirota's method. The inverse scattering [12, 13] and Hirota's 
method [10] results on N-dark-dark sohtons in the defocusing 2-CNLS model have presented 
only the degenerate case. 

Second, as in the AKNSi solution (5.25), one can show that the 1-dark-dark solitons 
(5.35) allow the imposition of the conditions (2.8) for the pair of fields {^f , ^^), i = 1,2, 
in order to satisfy the 2-CNLS equation (2.9). In fact, consider the parameters y'^ — 
and a complexification procedure as in (2.8) provided with a convenient set of complex 
parameters in (5.35), k'^ -> ik'', pf -> ipf such that y'' = 6^*"^% z'^ = ^\J^^jPjP] ~ i^'^Y^ 
4>q being real, so the equation (2.8) can be satisfied if 

(p+)* = //a,p-. (5.37) 

Notice that, as it was shown in the AKNSi case and its 1-dark sohton (5.26) particular 

reduction, in order to get 2-CNLS 1-dark solitons one must consider q = 1 in the both com- 
ponents {^j' , ^j)- Therefore, the relationship between the parameters pf determine clearly 
if these solutions will correspond to the defocusing {aj = —1) or to the mixed nonlinearity 
((7i = —(72 = ±1) 2-CNLS system, respectively. 

Third, let us discuss the extension of our results to the case of free field NVBC (3.29). 
In order to compute the relevant tau functions analog to the ones in (3.26)-(3.28) one must 
consider the expression (3.35) instead of (3.14). So, the tau functions become 

f+^(x,i) = (A<,|£;«e^_(x,^)*W/^[*W]-MA„)e-^ (5.38) 

= t) - Pt)f^ix, t) + (A,| [e-^''E%e^^]¥'^h[^^%' |A„) (5.39) 

fr\x,t) = (Ao|4l^e^(x,t)^(°)/i[^(°)]-i|Ao) e"^, (5.40) 

= {-p-{x,t)+p-)f^{x,t) + (A,| [e->^°4;^e>^°]^W/i[^(°)]-^ |A.) (5.41) 

f^{x,t) = {Xo\¥'^h[¥'Y'\Xo). , (5.42) 

where the group element Q^^{x,t) is defined in (3.31)-(3.32) and x(a;,t) is an ordinary func- 
tion. We have used ©'^ from cq. (3.32) and the decomposition (3.33) for Xo, as well as 
the properties (B.15)-(B.16) of the highest weight representation. The tau functions (5.38)- 
(5.42) exhibit the modifications to be done on the previous eqs. (3.26)-(3.28) in order to 
satisfy the NVBC (3.29). These amount to introduce the group element ©?. and the factor 

e-^. Then from (5.38)-(5.42) and the analog to (3.25) -^f ^ pt + ^ and *~ = p- - ^ 
one can get 



= p+(x,t) + ^; ^r=p-(a;,t)-^ (5.43) 

To To 
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f± = (A„| [e->^°£;i'j^e^°]*(°)/i[*(°)]-^ |A<,) ; = fo = (A„| ¥''^h[^^%^ |A<,) (5.44) 



In fact, the relations (5.43) can be formally obtained from the ones in (3.25) by making 
the changes pf — >■ pf{x,t), rf" — >■ r^, Tq — )■ tq. The soliton type solutions can be obtained 
following similar steps as the previous ones; i.e. for 1-soliton choose h'^ = e"-'''^''^^^ and the 
adjoint eigenvectors of £1^2 become precisely the vertex operator as in (5.28)-(5.29), then 
the expressions {k'^x + wH) of the tau functions (5.32)-(5.33) do not change. 

Next, let us describe the general properties of the solutions (5.43)-(5.44). i) Notice that 
the form of fo remains the same as the one in (5.32). ii) The modified DT and the related 
tau functions f± (5.44) will introduce some new parameters into the tau functions rf" in 
(5.33). In (5.33) the factor pf must be replaced by pf{x,t) and s^*^^) in general will depend 
on the index 'r and the space-time, i.e. s^^"^^ — )■ s^*^^'*''(a;, t). The explicit form one gets from 

W^'^IAo >. In the case of 1-dark-dark-soliton 

5^ ^'^^ , xf ^'^^ being constant 



s^(i'^\x,t) = {pf{x,t))-^ai < \, 

it is required an element Xo such that s^^'^''^\x,t) = ^'^'^ xf ^'^\ 
parameters. So, under the above modified DT one has the general form of the 1-dark-dark 
soliton 



(1 + (y])±^) + (-1 + (y|)±^) tanh{(A;«x + wH + log c^^VS}] , (5.45) 



where j = 1,2; = 



—a,- , b'l — —6,- and 



(5.46) 



The solutions (5.45) are the general 1-dark-dark solitons of the AKNS2 model and the two 
versions q = 1,2, to our knowledge, have never been reported yet. A complexified version of 
this soliton for q = 1, which is a solution of the 2— CNLS model, has recently been reported 
in [4, 14]. 

2-dcirk-dcirk solitons. In order to get the two-dark-dark soliton (A^ = 2) solution one 
must choose the group element — e"'^^''^''^'' e^'^^'^^'^'^K Then, following similar steps one has 



To 



= 1 + ci e'^i^+^i* + C2 e'=2^+"'2* + ciC2rfoe'='''+"'^*e'=^^+"'2*, 



+ s, e 



1,2 



(5.47) 
(5.48) 



where 



do = 




Si So 



Si S2 



j (^n) 



Si S2 
2 



Si s 



1 '^2 



SV S 



1 ■'2 



(Sr-S + )(S2 

2 + -"i 
k=l PkPk) . 



-do; (5.49) 



^kl (n = l,2)(5.50) 



The above process can be extended for A^— dark-dark soliton solutions, in that case the 
group element in (3.26)-(3.28) must take the form hi = e^i^i ('=i)e"2i^2('=2)___gajvVF« (fcjv) _ xhe 
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relevant tau functions follow 




NxN 



i = l,2. 



m,n = 1,2,3...A^ 



(5.52) 



(5.51) 



where 



NxN 



stands for determinant, c. 




and Wn, k. 



"n 



are the same as in (5.50). Notice that c^n = Cnm, Cnn = 0. We have omitted the index 
(q) in all the parameters above. Let us mention that we have verified these solutions up 
to A'^ = 3 using the MATHEMATICA program. The above solutions associated to the two 
(5.47)-(5.48) or higher order (5.51)-(5.52) dark-dark-solitons are in general non-degenerate. 
Notice that in the CNLS case the two— and higher— dark- dark solitons derived in [10] are 
actually degenerate and reducible to scalar NLS dark solitons. So, regarding this property 
our solutions resemble to the ones recently obtained in [4, 14] for the r-CNLS system. 

Let us remark that the free field NVBC (3.29) requires the introduction of more parame- 
ters into the N-dark-dark tau functions (5.51)-(5.52) through the modified DT. So, following 
similar steps to get (5.45) one has: i) the parameters entering the tau function Tq in 
(5.51) remain the same, ii) in (5.52) the factor pf must be changed to pfe"'j *. iii) the 
parameters entering the tau functions rf" in (5.52) must be changed as — )■ s^xf- 

So far, to our knowledge, general multi-dark-dark solitons in the AKNS model incorpo- 
rating the cases q — 1,2 have never been reported yet. In [4] the general non-degenerate 
N-dark-dark solitons in the r-CNLS model with mixed nonlinearity has been reported in the 
context of the KP-hierarchy reduction approach, and in [14] the r-CNLS model has been 
addressed in the algebro-geometric approach. 

5.3 Dark-dark soliton bound states 

So far, reports on mult i- dark- dark soliton bound states in integrable systems, to our know- 
ledge, arc very limited. Recently, it has been shown that in the mixed-nonlinearity case of 
the 2-CNLS system, two dark-dark solitons can form a stationary bound state [4]. Then, in 
order to have multi— dark— dark soliton bound states in the N-soliton solution (5.51)-(5.52) 
the constituent solitons should have the same velocity, i.e. denoting y„ = (j/„ ^ [0, 1]) 
then = ^f^kn = v (assume v > 0) for certain soliton parameters labeled by n = 1, 2, .... 
We show that the signs of the sum J2j Pj P] determine the existence of these bound states, 
for the positive sign bound states can be formed, whereas for the negative sign they do not 
exist. 

First, consider J2iPtP7 > 0) then from (5.50) it follows that A;„ = ^•'^"i— i: (^^ pt PiY^"^ fo'^ 
7/„ < . Therefore, in order for multi— dark— dark soliton bound states to exist the equation 



(|yn| - 1)^ 



V 



= c> 



(5.53) 



{\yn\ + l)^|\y^\ 



i Pi Pi 
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must give at least two distinct positive solutions for In fact, one has the following two 
solutions 



= l + - + c^^± , (5.54) 

where A = These exaust all the possible solutions with the condition (5.53). 

The right hand side of (5.54) with either + or — signs provides \yi,2\ > for any c > 0. 
The other possible solutions (/» ^ do not satisfy the above requirements, since they possses 
a term of the form \/12c^ + c'^ — A, which is imaginary for any positive value c. Therefore, 
two— dark— dark— soliton bound states exist in the s/(3)— AKNS system, and three— and 
higher— dark— dark— soliton bound states can not exist. These results hold for any value of 
the index q. Notice that in order to reduce to the 2-CNLS system the parameter relationship 

(5.37) must be satisfied. So, from (5.37) one has J2iPtpi — ~[^^^ + ^'^^] < (remember 
that under pf — >■ ipf the sum ^ipfp^ reverses sign); so, it is possible if ai — — (T2 (e.g. 
(Ti = —(72 = —1 for \pf\ > Iptl), this case corresponds to the 2-CNLS with mixed focusing 
and defocusing nonlincaritics. Thus, the two-dark-dark soliton bound state solution we 
have obtained here corresponds to the Manakov model with mixed nonlinearity. The case 
CTj = — 1 (i = 1, 2) defines the defocusing Manakov model which does not support multi-dar- 
dark-sohton bound states [4]. 

Second, the condition J2i pt Pi < will provide only a single positive solution for > 
in the equation ^ = t; > 0. So, one can not obtain two or more solitons with the same 
velocity and therefore bound states in this case are not possible. 



6 Mixed boundary conditions and dark-bright solitons 

One may ask about the mixed boundary conditions for the system (2.5)-(2.6), i.e. NVBC for 
one of the field components of the system, say , and VBC for the other field component, 
'^2- So, we will deal with the mixed boundary condition (3.4). Moreover, taking into account 
sequence of conditions as in (2.8) 

t^-it, [^tV = -P (^i ^7 = -P (^i A, (6.1) 

the system (2.5)-(2.6), for {/3i.p) being real, can be written as 

idtil;k + dlil;k + 2p\J2(^j\^j\^-7:{/3k.p)Uk = 0, k = l,2. (6.2) 

Precisely, the system (6.2) in the defocusing case (^oy = = ~1) and possessing the first 
trivial solution (3.4), i.e. p\ 7^ 0; = 0, has been considered in [26] in order to investigate 
dark-bright solitons which describe an inhomogeneous two-species Bosc- Einstein condensate. 
The system (2.5)-(2.6) with the mixed trivial solution (3.4) can be written as 

dt^t = ±dl^tT2\Y^^^i^- -pXpl\^t. (6.3) 



2 



d^^t = ±9^*rT2[E*;*7-^p;"pr-^^^2]^^ (6.4) 
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From the previous sections we can make the following observation: the vacuum connec- 
tions relevant to each type of solutions must exhibit the fact that dark solitons are closely 
related to NVBC, whereas bright solitons to the relevant VBC one. Since Q2 and Qi are some 
parameters satisfying 2Qi + 0,2 — '^PiPi, for simplicity we will assume Q2 = 0, D,i — PiPi 
in (3.6)-(3.7). Therefore, the connections (3.6)-(3.7) for the mixed (constant-zero) boundary 
condition (3.4) take the form 

^ sl^=E^'^ + ptEj^;>+p^E% (6.5) 
S''- = sl = E^'^+ptEj^^^+p^E% (6.6) 



The relevant group element *^^c given by 

^^l^e^''^.^'^'^, ■ (6.7) 



In order to construct the sohton solutions we must look for the common eigenstates of 
the adjoint action of the vacuum connections (6.5)-(6.6). So, one has 



[4,r±^,] = ±A;±r±0„ (6.8) 
sl,T±p,] = ±«;±r±^„ w^ = {k^f-ptp^ (6.9) 

where the vertex operators r±^2(A;='=) are defined in (C.3). We expect that these vertex 
operators will be associated to the bright-dark soliton solutions of the model. 
Let us write the following expressions 



[x£j^+i4,r^J = <p+{x,t)r0,; <p+{x,t)^k+x + wn; (6.10) 
ce^j , r_/32 = -ip~{x,t)T_is2; (p~{x,t) = k~x + w~t. (6.11) 

The 1-dark- bright sohton is constructed taking h in (3.26)- (3.28) as 

^ g7+r+^2 e^-r_^2 ^ 7± = constants, (6.12) 

where the vertex operators in (C.3) have been considered. Notice that due to the nilpotency 
property of the vertex operators, as presented in the appendix C, the exponential series must 
truncate. So, replacing the group element (6.12) in (3.26)-(3.28) one has the following tau 
functions ^ 

f± = e'^^-^'j+r <\o\E^^It+^,T_^,\Xo> (6.13) 
= e±'^*7± < A<,|4'^^r±^,|A„> (6.14) 
To = l + e^^-^"7+7- < A<,|r+;3.r_^.|Ao> (6.15) 



^There exist other eigenstates [e^^ , V^] = XV^; [ij^ , V^] = (-1)«-U(A2 - pl^^^V^^, wherepg = 
4:PiPi and Vg^ {q = 1, 2) is given in (C.5). However, these eigenstates are related to purely dark-solitons for 
the first component e.g. if one takes the group element h = e -^i , it will not excite the second component 
since the matrix elements of type < Xo\E^^^^V^^ \Xo > vanish. 
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These tau functions resemble the ones in (5.32)-(5.33) for the 1-dark sohton (5.35) and the 
eqs. (4.8)-(4.9) for the 1-bright sohton (4.10) components, respectively. The matrix elements 
in (6.13)-(6.15) can be computed, and will depend only on the parameters fc^, pf. So, one has 
six independent parameters 7^, /c^, pf associated to the these tau functions . Let us write the 
tau functions in the form fo = 1 + cq e 
^± g±fe±a;±«;±t^ where = {k^y - Pi Pi', Co 



k~^x+w^t„—k x—w t ^± 



„it „A;+x-+ui+t „— fc x—w t ^it 
(X o c 5 ' 2 

k+ 



a'^Pi -a Pi 

l){k'^k~ + Pi Pi). So, using these tau functions in (3.25) one has 



ai ^ p^;b+b- = (^)(F 



Pi 



± 



pf ± ( — ) 1 + tanh 



{k'^ — k )x + (u!"*" — w )t 



± 



2v^ 



(6.16) 
(6.17) 



This is the 1-dark-bright soliton of the sl{3) AKNS model (6.3)-(6.4) (for Q2 = 0). Notice 
that this solution has six independent complex parameters, say a~,b~,k^,pf. 

The construction of the 2-dark-bright solitons follows similar steps. The group clement 
h = e'^'i ^+i^2(''i ) e'^'i ^-(^2^^! )e^^ ^+P2'^^2 ) e^"^ ^-P2^^2 ) does the job. We record the relevant tau 
functions 



To = 



ft 



1 + 

2 

m,n=l 
2 



m,n=l 
2 



,± „±kmX±Wr„.t 



where = {k^y - pfpi ; c 



m=l 

±\2 



m=l 



Pi 



Co 



d+d^ d^ 



(6.19) 
(6.20) 



^1 ^2 Pi 



"12(^1-^2 ^2 +P1P1 ) _|_ '^22(^2-^2 )(fe2 ^2 +P1P1 ) 

a;"i(fe;*"-fcf )(fc+fc;"+p+pj")) a2-^(fcf-fc+)(fc;"fc++p+p;") 



0; bTh 



(k+-k-){k+k-+p+p-) , 



"l2"2l(fel -^2 ) (^1 -^2 ) (fcl ^2 -fel ^2 ^2 +Pl Pi ^2 +Pl Pi ) ■ h_ 



(k+-k^f{k^-k+)[k+-k:^){k+-k^)'^pl{k+k^+p+p^){k+k^+p+pl)'bl 



"12^1 (fcl -^2 )(*^1 ^2 +Pl Pi ) ■ 

a^jA;^ )(fc^fef +p^pf )' 



^(^1 ^2 )«: 



± 



W (pfF {kt-k7)ikf-kf_^) ' 

These tau functions resemble to the ones in [24] provided for the 2-CNLS model. The 
generalization to N-dark-bright solitons require the group element h to be 



fi — Qyt'^+02(kt) e^i ^-132(^1 )g7jr+^2('=J) e^2 r-^2('=2 )....e'^^r+;32(A;+) g7jyr_^2('=jv). 



(6.21) 



7 Discussion 

We have considered soliton type sohitions of the AKNS model supported by the various 
boundary conditions (3.2)-(3.5): vanishing, (constant) non-vanishing and mixed vanishing- 
nonvanishing boundary conditions related to bright, dark and bright-dark soliton solutions. 
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respectively, by applying the DT approach as presented in [15]. The set of solutions of the 
AKNS^ system (2.5)-(2.6) is much larger than the solutions of the r-CNLS system (2.9). A 
subset of solutions of the AKNS^ system, (2.5)-(2.6) for r = 2 and (5.3)-(5.4) for r = 1 , 
respectively, solve the scalar NLS (5.6) and 2— CNLS system (2.9), under relevant complex- 
ifications. 

Moreover, the free field boundary condition (3.29) for dark solitons is considered in the 
context of a modified DT approach associated to the dressing group [29], and the general 
N-dark-dark soliton solutions of the AKNS2 system have been derived. These sohton compo- 
nents are not proportional to each other and thus they do not reduce to the AKNSi sohtons, 
in this sense they are not degenerate. Wc showed that these solitons under convenient com- 
plexifications reduce to the general N-dark-dark solitons derived previously in the literature 
for the CNLS model [4, 14]. In addition, we have shown that two— dark— dark— soliton bound 
states exist in the sZ(3)— AKNS system, and three— and higher— dark— dark— soliton bound 
states can not exist. These results hold for any value of the index q. In the case of reduced 
2— CNLS when focusing and defocusing nonlinearities are mixed, this result corresponds to 
2-dark-dark soliton stationary bound state [4]. 

In the mixed constant boundary conditions we derived the dark-bright solitons of the 
sZ(3) AKNS model. These solitons under the complexification (6.1) reduce to the solitons 
of the 2-CNLS model (6.2) which will be useful in order to investigate dark-bright sohtons 
appearing in an inhomogeneous two-species Bose- Einstein condensate [26]. 

Another point we should highlight relies upon the possible relevance of the CNLS tau 
functions to its higher-order generalizations. We expect that the tau functions of the higher- 
order CNLS generalization are related somehow to the basic tau functions of the usual 
CNLS equations. This fact is observed for example in the case of the coupled scalar NLS-I- 
derivative-NLS system in which the coupled system possesses a composed tau function de- 
pending on the basic scalar NLS tau functions [40]. 
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A sl{2) matrix elements 

The commutation relations for the sl{2) affine KacMoody algebra elements are 



±2E. 



(m+n) 



E. 



(m) ^(n) 



(A.l) 
(A.2) 

(A.3) 
(A.4) 



The central extension ensures highest weight representations (h.w.r.) of the affine algebra 
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(see e.g. [15]). So, in the h.w.r. {|Ao >, |Ai >} one has the following relationships 

£^f|A„> = (A.5) 

Ei™^|Aa> = 0, m>0 (A.6) 

i/M|A^> ^0, m>0; (A.7) 

/^^°^|A„> = 5a,l\\a>, (A.8) 

C|A„> = |A„> (A.9) 

where a = 0,1. The adjoint relations = i?:^"™-*, )t = _f/'(-"») allow one to know 

their actions on the < Xa\. Next, consider the vertex operators 



^ ^^"^T + e, (^^^)^/^<oC^}; g = 1, 2; (A.IO) 



7 + (72-p2)i/2-^- 



where = (—1)^ ^. 

The following matrix elements can be computed using the properties (A.5)-(A.9) 

(A„|l-^|A„), = e,^^^^ (A.ll) 

7 



(A„|4"V«|A„) = ^—^-^^(-r^-p^'. (A.12) 

The matrix element (Ao| V^^(7i, p)V''^(72, p) |Ao) can be computed by developing the products 
and keeping only non-trivial terms, then one makes use of the commutation rules to change 
the order, and eventually to get some central terms C. The double sum can be simplified to a 
single sum and each term can be substituted by power series like J2^=q x"' — (iz^' Sn^i ~ 
(iz^, and E~=i^a;" = J]^- So, one can get 

<A„|V'(7.,P)V"(72,P)|A„) = 12 + 2A-(7„72)lp^||^ + 

[ ^'^"^^'"' + 1]; ,= 1,2 (A.13) 

7i72 

where A-.(,..,.) . , 
\/7? - P' 



In order to prove the nilpotency property of the vertex operator V'^{ji, p), when evaluated 
within the state |Ao >, it is convenient to write (A.13) in the following Laurent series 
expansion 

(A.|V-(,,.p)V«(7..p)|A„) = -6/y^v/^( ^^^ f^ )^(^)' X 

72 - P 71 + 72 
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From this, it is clear that 



IO72 (71 - 72) , 15(67! + P^) (71 - 12)' 



+ 



4! (7i-p2) 5! {ll-p'Y 
42072(27! + (71 -72)=^ 
ill - P'f 7! 



6! 



+ 



Um (A,| y^(7i, p)T^n72, p) \K) ^0, g = 1, 2. 

71— >72 



(A.15) 



(A.16) 



B The affine Kac- Moody algebra sl^{C) 

In the following we provide some results about the affine Kac-Moody algebra Q = 5/3(6*) 
relevant to our discussions above. We follow closely [41, 18]. The elements of the 5/3(6*) 
Lie algebra are all 3 x 3 complex matrices with zero trace. Consider the corresponding root 
system A = {±ai, ±0:2, ±03}, such that the three positive roots are ai, i = 1,2, 3, with «„, 
a — 1,2, being the simple roots and eta = cti + 0:2. We choose a standard basis for the Cartan 
subalgebra H such that 



( 1 

















-1 






1 


i) 


vo 







VO 








(B.l) 



and the generators of the root subspaces corresponding to the positive roots are chosen as 



E. 



+ai 




+012 



/ 


V / 



n 



E. 



+0,3 




(B.2) 



For negative roots one has E_a — (E+a)'^- The invariant bilinear form on s/3(C), {x\y) = 
ti^xy); x,y G sls{C) induces a nondegenerate bilinear form on l-L* which we also denote by 
(• I •). This definition allows one to write 



(ailcKi) = 2, (q;2|q;2) = 2, {ai\a2) — —1 
On the other hand, the generators T^"^^ = {h["^\ h'^\ Ej^^}, where m G 2 and 



together with the central C and the 'derivation' operator D ([L',T("*)] 
basis for sl3{C). These generators satisfy the commutation relations 



m 



pirn) pin) 

pim) pin) 
p(fn) Tp{n) 



= m {aa\ab) C 6m+n,0, 
^ pim+n) 



(B.3) 

a e A, 
form a 



(B.4) 
(B.5) 
(B.6) 
(B.7) 
(B.8) 



24 



pim) jpin) 
in) 

-C12 



, E. 



= 0, 



(B.9) 

(B.IO) 
(B.ll) 



where a,b = 1,2, i = 1,2,3 and m, n e Z. The remaining non- vanishing commutation 



An) 



rpi-m) rni-n) 



relations are obtained by using the relation E^™) , 

In this paper we use the homogeneous Z-gradation of slziC) which is defined by the 
grading operator such that 



shiC) = ^ Qm, [Gm, Gn] C Gm+n, 
meZ 

Where Gm = {x E sl^^C) \[D , x] = mx; m E Z}. 
The subspace Go is a subalgebra of 5/3(6*) given by 

Go^CHi®CH2® CE^^^ ®CC®CD 

and for the subspaces Gm ijn 7^ 0) we have 

Gm = CH^"^ ®CH^^^ (BCE'^^^ ®CE^^^ (BCEt'^ (BCE^:::'! 



(B.12) 



as 



ai 



CE. 



(m) 
—0.2 



(B.13) 



(B.14) 



We use in the paper the fundamental highest weight representation | Aq ) , satisfying 



Ao) 



(B.15) 



i/W|Ao)=0, E(°)|Ao)=0, C\Xo) 

ior a,b — 1, 2, and a E A. Such state is annihilated by all positive grade subspaces 

g^\Xo)^0, m>0, (B.16) 

and all the representation space is spanned by the states obtained by acting on |Ao > with 
negative grade generators. This representation space can be supplied with a scalar product 
such that one has 



C7t = C, 



= D. 



rp—m 
E-a : 



It follows from (B.13) and (B.14) that 



(^m)^ — G—m-i 



and therefore 



< AqI G-r 



0, 



m > 0. 



(B.17) 
(B.18) 



(B.19) 



(B.20) 



In addition to the subalgebra go it is also convenient to consider two additional subalge- 



bras 



G<Q — ® Gm^ 

m>0 



^>0 — ^ G-m- 

m>0 



(B.21) 
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These subalgebras and the corresponding Lie groups play important role in the DT 
method. 

The next relationships are useful in the AKNS^ (r = 2) model construction. The special 
element E^^^ in the basis presented above can be written as 

£;(0 ^ i(ii-(0 + 2Hi^^), [D , = lE^^K (B.22) 
Then the matrix E^^^ becomes 




(B.23) 



The roots entering in the AKNS2 construction are 

/3i = as = «i + "2; f^2 
Moreover, the following commutation relations hold 



a2. 



Z+m,05 



E. 



Ml 

(m) 



±E. 



±/3iT/32 



(l+m) 

0, 



77(") 



(m+n) 
±/3i ' 



I 77i("i+n) 



a = 1,2, 
= 1,2. 



Trim) rp(n) 
-"1 ' ^±132 

Trim) rp(n) 



77.(") 



m+n) 



7T.(m+n) 
+^±/3i=F/32 



p(m+?i) 



= ±2K 



(m+n) 

'±/32 



(B.24) 

(B.25) 
(B.26) 
(B.27) 
(B.28) 
(B.29) 
(B.30) 
(B.31) 



C s/(3) matrix elements 

Consider the vertex operators associated to bright sohton solutions 

+00 +00 
Fj- E E P"4r^; i = l'2; uj,pje€. (C.l) 

n=— 00 n=— 00 

It can be shown that they are nilpotent, i.e. Fj — 0,G'j — 0. The matrix element 
{Xo\ FjGk\Xo) can be computed by developing the products and keeping only non-trivial 
terms, then one makes use of the commutation rules to get the central term C. The double 
sum can be simplified to a single sum, which provide the power series Z^^i nx" = (X^p-- 
So, one has 

(Ao| FjGk \Xo) = - — ^-^-^-2^j,k (C.2) 
W - Pk) 
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Let us consider the deformation of the vertex operators F2, G2 as 

+00 ^ 

-i2 



r±,Ak^)= E (^)-"[^^4?-p?4tm], w^^{k^r-ptp^. (C.3) 



It is a direct computation to show the nilpotency of these operators, i.e. T^.^^ = 0. 
Similar computations to the one in (C.2) provide the following matrix element 

(A„|r,.(*-)r_,.(r)|A.> = (C.4) 

Consider the vertex operator analog to the one in (A. 10) 

00 -1 + 



where = (—1)^^''^ and pg = Api pi . The next matrix element computation follows similar 
steps to the one performed to arrive at (A. 13), except that one must take into account the 
s/(3) commutation rules. So, one has 

{K\Vl(Kp«)Vl{\2,P,)\K) = i(P + 2-ff(Ai,A2)l^^j|^^^ + 

l^^^ + l]}; ,= 1.2, (C.6) 

A1A2 

where KoandK are given in (A. 14). Since this two-point function, except for an overall 
constant factor, is similar to the one in (A. 13) one can use the relationships (A.15)-(A.16) 
to show that the operator V^^(Ai,po) is nilpotcnt. 

The vertex operator generating the dark-dark soliton solution becomes 

00 2 2 

W^iKpt,,) = E {(k' - 4EP.^Pr)-"/^[e,]"hi^i"^ + s,Ht^ + E-4)4:^ + 

n=— 00 i=l i=l 

(C.7) 

± ^ Tpf ^1 ^ 1 pjpi ± ^ 1 P1P2 

kTe,ik^-4. Ej Ptpjy/' ' 2 ' 2 Ptpi ' 2 E^ pIp7 ' 



(A^ - pI) 



Notice that the vertex operator (C.7) reduces to the one in (C.5) in the limit 0. 
The nilpotent property of this vertex operator can be verified as follows 

(A,| W\h, pt2)W\h, pt,) \\o) = (^)^^ (^1 + - X,) + ..\ (C.8) 



where Xi = ^ kf — AS, X2 = y k2 — AS, S — Y,j Pj^Pj ■ the limit X2 — >■ xi (or k2 ^ ki) the 
r. h. s. of eq. (C.8) vanishes. 
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